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Abstract: Using a numerical model we study the frequency de- 
pendence of mode coupling gain due to stimulated thermal Rayleigh 
scattering in step index, Yb doped, fiber amplifiers. The frequency at 
the gain peak is shown to vary with core size, doping size, population 
saturation, thermal lensing, fiber coiling, direction of pumping, 
photodarkening, and pump noise spectra. The predicted frequencies 
are compared with measured values whenever possible. 
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1. Introduction 



In earlier papers we described the physical process that causes modal instability [1] in 
high power fiber amplifiers. A key aspect of this stimulated thermal Rayleigh scattering 
(STRS) process is a frequency offset between the two coupled modes, usually LPqi 
and LPn. We have incorporated the physics of this STRS process in a detailed, highly 
numerical model [2] that can compute mode coupling gain as a function of the frequency 
offset. A typical gain curve is shown in Fig. 1. In [1] we argued that the frequency of 
peak gain could be roughly estimated as the inverse of the thermal diffusion time across 
the core radius. This implies the peak gain for typical large mode area fibers should lie 
in the frequency range 200-5000 Hz. 
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Fig. 1. Mode coupling gain plus lasing gain versus frequency offset for a large 
mode area fiber amplifier. The coupling is from LPoi to either LPn or LPo2- The 
peak value of gain and its position vary with amplifier parameters. 

Semi analytic versions of our model of STRS have been presented by Hansenei al. 
[3,4]. Hansen et al. also computed gain versus frequency offset, but with the assumptions 
that the transverse heat profile is proportional to the transverse irradiance profile and 
that the mode profiles are constant along the fiber. These simplifying assumptions 
permit the development of a compact model, but such models miss some of the subtleties 
of real amplifiers. For instance, in real amplifiers the shape of the mode coupling gain 
g(z) is altered by any change in the transverse heat deposition profile along the fiber, 
and depends not only on the core refractive index step and the doping profile, as found 
by Hansen et al., but also on the bend radius, linear absorption of the signal due 
to photodarkening or other processes, thermal lensing, and with total thermal load if 
cooling is asymmetric. However, perhaps the most important difference between the 
simplified models and our more detailed model is related to changes in the transverse 
profile of the population inversion, or transverse hole burning, which dramatically alters 
the shape of the oscillating heat source responsible for STRS as the balance between 
pump and signal powers changes. This is most noticeable in co-pumped amplifiers where 
the pump power falls while the signal power rises as they propagate down the fiber. 
However, hole burning is also important in counter-pumped amplifiers, even though 
the ratio of pump to signal powers is more nearly constant along the fiber. Population 
depletion affects the magnitude of the mode coupling gain as well as its frequency profile, 



but its influence on the magnitude of gain will be the topic of other papers. 

In this report we present modeling results for amplifiers operated at the gain peak 
and near the instability threshold. We expect our steady-periodic model to work well 
near threshold although it may be less useful for operation well above threshold. Here 
we are primarily concerned with how the frequency that gives maximum gain near the 
instability threshold (Fm) changes with fiber design and operating conditions. Besides 
being an important test of our model, this information will be useful in future modeling 
since it reduces the scope of the two dimensional search required to find Fm at the 
threshold. We also anticipate that measured values of Fm will be useful in diagnosing 
causes of anomalously performance such as low instability thresholds. As we will show, 
different problems can have different influences on the frequency. Even so, it is unlikely 
that in practice frequency information alone will uniquely identify such problems. 

2. Frequency vs mode field diameter (or effective area) for heat profile 
proportional to irradiance profile 

In this section we show how Fm varies with core size, using the simplifying assumption 
that the heat profile mimics the irradiance profile. This is the same approximation 
used by Hansen et al. [3]. It gives a preliminary estimate of Fm that we will improve 
on in the next section. For this preliminary computation we use a small linear signal 
absorption inside the doped region, with no pump absorption. This produces a heat 
profile proportional to the irradiance profile in the doped region. All the power lost to 
linear absorption is assumed to appear as heat that causes mode coupling via STRS. 
We use a step index fiber with numerical aperture (NA) of 0.054 and a 976 nm pump. 
Our peak frequencies for the 20, 40, and 80 flm diameters agree closely with those of 
Hansen et al. [3] even though we use a 1040 nm signal while Hansen uses 1030 nm. This 
is expected because Fm depends only on the shapes of the two modes and they do not 
change significantly over this wavelength range. Our results are summarized in Fig. 2. 
The frequencies are well approximated by the equation 

F M = 1.71xl0 6 /A eff . (1) 

where A e g is the effective area of the LPoi mode in square microns and the signal 
wavelength is 1040 nm. This result can be compared with a frequency estimate based 
on the thermal diffusion time for a heated cylinder of radius r e g- 

F=^~ (2) 

which for silica thermal parameters (K = 1.38 W/m-K, C = 703 J/kg-K, p = 2201 kg/m 3 ) 
gives 

F = 8.9xl0 5 A e 2 ff . (3) 
The two frequencies computed using Eqs. (1) and (3) are equal if 

A eff =1.9r e 2 ff . (4) 

Defining the mode field radius r m f by 

A e & = 7tri f (5) 

means 

r e ff= 1.28/mf, (6) 



a reasonable result. 

In the computations above we fixed the numerical aperture at NA = 0.054. Chang- 
ing the NA changes the shape of the modes and thus the frequency. Hansen et al. [3] 
showed the frequency increases as the V parameter, or equivalcntly, the NA increases. 
An increase in NA tends to compress the modes so it is no surprise Fm increases slightly. 

3. Realistic heat profiles 

In real amplifiers the heat profile does not match the product Eq\E\\ as assumed above 
because the upper state population is depleted by varying amounts across the core. The 
degree of depletion varies with the ratio of signal to pump irradiances, so it also varies 
along the length of the amplifier. Figures 3 and 4 show how the signal and pump powers 
vary along the fiber for high efficiency co- and counter-pumped amplifiers. For the co- 
pumped amplifier the pump and signal powers are approximately equal midway along 
the fiber. For the counter-pumped amplifier they are approximately equal along most 
of the fiber. Near the input of the co-pumped fiber the pump is much stronger than the 
signal so there is little population depiction, or hole burning, and the approximation of 
the previous section is nearly met, and Fm is given by Eq. (1). The time-averaged heat 
profile at this location is shown in blue in Fig. 5, while its oscillatory component is shown 
in blue in Fig. 6. The corresponding heat profiles midway along the fiber are shown as 
green curves in the two plots, and the profiles near the output are shown as red curves. 
For the same amplifier with counter-pumping the profiles are nearly constant along the 
fiber and are approximately the same as those near the midpoint of the co-pumped 
fiber. 

Clearly, transverse hole burning strongly alters the profile of the anti symmetric, 
oscillatory part of the heat deposition. The assumption of the previous section which 
ignored hole burning as unimportant is highly dubious for real amplifiers, whether co- 
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Fig. 2. Frequency of peak gain (Fm) vsAgg of mode LPoi for a constant numerical 
aperture of 0.054 and a 1040 nm signal. The heat profile matches the irradiance 
profile over the doped region in this case. The circles are found using our numerical 
model for core diameters of {25 35 40 50 60 70 80 90} flm with corresponding 
effective areas of {380 645 805 1180 1630 2155 2755 3420}, while the solid line is 
given by the best fit inset equation. 
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Fig. 4. Signal and pump power distributions along typical counter-pumped am- 
plifier. 
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Fig. 5. Transverse cuts through the normalized profiles of the time averaged heat 
at three z locations along the 1.2 m long, co-pumped fiber. The fiber has rf core = 
d(j ope = 35 /lm and is designed for high conversion efficiency with iVyb = 3 x 10 25 
m~ 3 and d ciad = 140 /lm. It is co-pumped with 425 W, seeded with 20 W of 1040 
nm signal in LPqi and 0.2 W in LPn- The powers evolve as shown in Fig. 3 
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Fig. 6. Transverse cuts through the normalized profiles of the anti symmetric, 
oscillatory portion of the heat at three locations along the co-pumped fiber. The 
amplifier is otherwise the same as in Fig. 5. 



or counter-pumped. The outward displacement of the oscillatory heating means the 
asymmetry persists longer before it decays away by thermal diffusion, and this means 
Fm is reduced as hole burning strengthens. For the co-pumped fiber most of the mode 
coupling gain occurs near the mid point, so the heat distribution at that location largely 
determines Fm- The value of Fm for the counter-pumped fiber is approximately the same 
because the heat profile over the full fiber length is nearly the same as at the midpoint 
of the co-pumped fiber. 

Full model runs using computed population inversions are used to predict realistic 
frequencies for the 20, 40, and 80 flm diameter, 1.2 m long, co-pumped amplifiers of 
Fig. 2. The pump cladding diameter is scaled with the core diameter to maintain similar 
pump to signal ratios and similar conversion efficiency. The results are shown as the 
symbols in Fig. 7. For the (20, 40, 80) jj,m diameter cores the simple model gives (6130, 
2125, 620) Hz while the full model gives (5350, 1750, 475) Hz, making the ratios (0.87, 
0.82, 0.77). 




Fig. 7. Comparison of Fm computed using the simplified (solid line) and full models 
(circles) for core diameters of 20, 40, and 80 flm. All three are 1.2 m long, with 
NA = 0.054, co-pumped at 976 nm, with X s = 1040 nm, and operating near the 
mode instability threshold. 



4. Dependence on doping diameter 

It can be anticipated that if the Yb 3+ dopant is confined to a smaller diameter than the 
core index step, the influence of hole burning on Fm will be less than for a fully doped 
core because heating will be confined to the doped zone and cannot broaden as much. 
Keeping other fiber properties constant and varying only the doping diameter, Hansen 
et al. [3] calculated the frequency for three doping diameters with a fixed refractive 
index profile. They showed Fm increased with decreasing doping diameter. However, 
they used heat profiles proportional to the irradiance profiles. We have recomputed the 
frequency response using our numerical model for one fiber design, varying the doping 
diameter and adjusting the pump power as necessary to achieve threshold operation. 
Figure 8 shows our results. In all three curves the power is near the mode instability 
threshold so the fraction of signal power in LPn is approximately 0.3% at the output 
for the frequency of highest gain. The fiber is the step index equivalent of LPF45 



(^core = 81 Jim, d clad = 255 fim, L = 1.2 m, n core = 1.45015, « c i a( j = 1.45, AA = 0.0209, and 
V = 10.2) with signal input in LPoi=10 W, signal input in LPn = 10~ 7 W. The gain 
peaks corresponding to doping diameters (63, 72, 81) /Xm lie at frequencies (500, 430, 
380) Hz. There is probably no deep significance to it, but the product of Fm and ifdope 
is nearly constant for these three cases. 

5. Dependence on V parameter 

A lower value of the numerical aperture would allow the modes to expand slightly into 
the cladding, and this should reduce Fm- This effect has been verified by Hansen et 
at [3]. under the assumption of no hole burning. 

6. Dependence on thermal lensing 

When the signal light is largely confined to the fundamental mode, as it is below and 
near threshold, the temperature peaks at the center of the core, and this causes thermal 
lensing, or self focusing of the signal light. The result is a constriction of the mode 
profiles which might be expected to cause an increase in Fm- The thermal lensing effect 
is stronger in larger diameter fibers. For example, it can be quite significant in an 80 
jllm diameter fiber operating at several hundred watts, but it is much less noticeable for 
a 20 fim diameter fiber. In Fig. 9 we show the effective area of the fundamental mode 
versus position along the fiber for two cases, a short and a long fiber with identical 80 
jum diameter cores and equal pump cladding sizes. The length of one is 1.2 m, the other 
is 2.4 m. The core is doped over its full diameter, but at a density of 3 x 10 25 m -3 for the 
short fiber and half that for the long fiber. This makes the hole burning similar for the 
two, but the heat per length is halved in the longer fiber, reducing its thermal lensing. 
As the figure shows, lensing reduces A e ff near the mid point by 13% for the shorter fiber 
and 7% for the longer fiber. 

If Fm scaled in proportion of 1 /A e ff this would cause a frequency shift of approximately 
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Fig. 8. Mode LPu fraction at fiber output versus frequency offset for varying 
doping diameters in co-pumped fiber. The core diameter is fixed at 81 microns 
and the doping diameters are 81, 72, and 63 |im; 1=1.2 m; rf c i a( j = 255 fim; A^b = 
3 x 10 25 m~ 3 . In each case the pump power is adjusted to the mode instability 
threshold. 



40 Hz between the two fibers. In fact, the modeled frequency shift is only 20 Hz, with 
the longer fiber having a lower frequency, as expected. The heat and LPn gain profiles 
of the 1.2 m fiber are shown in Fig. 10. The shift reduction is partly due to the fact that 
the gain peaks earlier in the fiber than the thermal lensing, and partly due to the fact 
that the gain is distributed over sections of fiber that have less lensing than the mid 
point. The result is that the effective frequency reduction is less than guessed based on 
the maximum reduction in A e g. Additionally, the smaller beam in the more strongly 
lensed short fiber leads to hole burning, and this also reduces Fm by small amount, 
bringing its shift closer to that of the long fiber. The lesson is that the influence of 
thermal lensing on Fm is less than might be guessed based only on the reduction in A e g-. 
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Fig. 9. A e ff for an amplifier with rf core = ddope = 80 flm. Other parameters: NA = 
0.054, rf clad = 376 jllm, LP i seed power is 10 W, LP n seed power is 10~ 16 W. The 
rapid oscillations are caused by launch of seed light with a profile equal to the 
low power mode instead of a self consistent thermally lensed mode. 



7. Dependence on coiling radius 

When a fiber is coiled all its modes are pushed toward the outside of the coil. One 
example is shown in Fig. 11. Not surprisingly, the values of Fm for the two orientations 
of LPn are shifted from the value in the same fiber without coiling, and the degeneracy 
between the two is broken. Figure 12 shows how Fm shifts with a changing coil radius. 
These frequencies were computed using the assumption that the heat profile matches 
the irradiance profile, as was done in Fig. 2. They are meant as a qualitative illustration 
of the increase in Fm with the stronger mode compression caused by tighter coiling. 

8. Dependence on pump wavelength 

If the pump light is detuned from the absorption peak at 976 nm, the reduced pump 
absorption coefficient will lead to lower upper state populations. This has the same 
effect as increasing the pump cladding diameter. It tends to increase hole burning and 
thus decrease Fm- 




9. Dependence on photo darkening 

The influence of photo darkening on Fm depends strongly on the photo darkening model. 
If we use a model with signal absorption that is uniform across the core, the frequency 
should be pulled toward the value from Eq. (1), that is toward higher frequency than 
the full model. In addition, the added thermal lensing caused by absorptive heating will 
also pull toward higher frequency. However, if we use a photodarkening model in which 
the absorption is larger in areas with higher upper state population [5] , the absorption 
will be strongest near the outer edge of the core, and this will tend to reduce Fm- We 
have not tested the latter model, but we have verified the shift to higher frequencies 
with the former. 

10. Comparisons with measurements 

A major problem encountered in comparing our computed values of Fm with measured 
frequencies is the general lack of information on the noise spectrum of pump or seed 
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Fig. 10. Blue curve is gain of LPn, including laser gain and mode coupling gain 
versus z for the same fiber as in Fig. 9. A e ff for an amplifier with rf core = Jdope = 80 
;Um. Red curve is the total heat deposited at location z in normalized units. 




Fig. 11. Irradiance profiles of LPoi (left image) and the two orientations of LPn 
in a step index fiber with numerical aperture 0.054, core diameter of 40 |im, coiled 
to a 100 mm radius. LPn (center image) is the mode with lobes out of the bend 
plane; LPn/ (right image) is the mode with lobes in the bend plane. 



modulation. Both noise spectra are likely to fall in strength with increasing frequency. 
This means the observed frequency would tend to shift to lower values. Unfortunately, 
to date the reported frequency measurements have not included measurements of either 
the pump amplitude or spectral noise or the signal amplitude noise. It seems likely that 
pump or seed modulation is present, however, because both produce signature effects 
that seem to be present in measured performance. Thus there is circumstantial evidence 
of pump or seed modulation. We hope that in future experiments the fiber inputs will 
be carefully measured and documented. 

Although data is as yet not well enough documented for rigorous comparisons, we 
can test whether our predicted frequencies and trends are in reasonable agreement with 
measurements. We discuss published examples below. 

10.1. Example #1 

Karow et al. [6] found F M « 1500 Hz in an NKT fiber model DC-400-38-PZ-Yb. Prop- 
erties of this fiber include d C0Te — c/dope = 38 fim, af c i a d = 400 fim, NA = 0.03, V = 3.36, 
A = 1064 nm, c/ mt =28 jttm (extrapolated from DC-200-40 because 400-38 is not listed on 
NKT products), L= 1.45 m, 9 W signal in, efficiency of conversion of absorbed pump 
of 0.80, a bend radius=20 cm, counter-pumped at 976 nm. They measured a signal 
noise spectrum that falls off with frequency roughly as 1 //, but it has multiple promi- 
nent spikes, making meaningful comparison with our model questionable. The measured 
threshold is at an absorbed pump power of 94 W which produces 77 W of output signal. 

We compute A e ff of 990 /im 2 for LPoi. The unusually low measured threshold power 
means there is insufficient heat load to cause significant thermal lensing. From the 
formula of Eq. (1) Fm would be 1730 Hz. The lower measured frequency of w 1500 Hz is 
consistent with the frequency reduction expected from the population saturation effect. 
However, a more quantitative comparison requires a more thorough knowledge of the 
input signal and pump noise spectra. 
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Fig. 12. Frequency of the gain peak for different coiling radii for coupling between 
LPoi and LPin (circles), and between LPoi and LPn (crosses). The fiber is the 
same as in Fig. 11. 



10.2. Example #2 

Otto et al. [7] measured a frequency of 350 Hz for a counter-pumped NKT fiber model 
LPF45. Parameters of this fiber are d com = 81 fim, c/^ope = 63 |Um, a! c i a( j = 255 jUm, L = 1 .2 
m, of m f = 67.65 /Xm, A e ff = 3595 fim 2 for LPqi. The measured threshold power was 200 
W for the output signal. 

We modeled a counter-pumped LPF45 simulated by a step index with V — 10.2 to 
match A e ff for LPoi , and found Fm = 490 Hz. Other model inputs included: 10 W signal 
seed in LPoi and 10~ 7 W frequency shifted light in LP11; 300 W pump input; 269 W 
signal output. 

A meaningful comparison between measurement and model is probably not possible 
in this case because the measured spectra are consistently highly structured, with three 
peaks at 300, 350, and 400 Hz, the most prominent being 350 Hz. We infer that the 
spectra do not reflect amplification of white noise but some other source of modulation, 
most likely either of the pump or the signal seed. The observed frequencies arc not 
inconsistent with the modeled value of 490 Hz. A meaningful comparison will require 
careful measurements of the input pump and seed spectra. 

10.3. Example #3 

Otto et al. [7] also measured the frequency for NKT fiber model DC-285/100 with 
o'coi-e = 100 jitm (step index), afdope = 100 fim, d c \ad = 285 jUm, d m f = 76 fim, A e g = 4536 
fim 2 , L = 1.2 m. The measured frequency was 220 Hz for a counter-pumped fiber at a 
threshold power of 100 W of signal output. The frequency according to Eq. (1) would 
be 375 Hz. 

We model a co-pumped, step index fiber with NA = 0.054, signal seed of 10 W in 
LPoi, 10~ 5 W of frequency shifted LPn, and 150 W of unmodulated pump. Our model 
predicts Fm = 325 Hz for a co-pumped version which is substantially higher than the 
measured 220 Hz. A small amount of photo darkening might reduce the frequency to 220 
Hz, or maybe the noise is 1/f like and so produces a lower frequency. Thermal lensing 
reduces A e ff by only 2% at this power level, so that should not alter the frequency 
noticeably. Again, careful characterization of the input light is needed. 

10.4- Example #4 

Ward et al. [8] measured a frequency of approximately 2000 Hz (with second, third, and 
fourth harmonics of 2000 Hz) for a fiber with d cole = 39.5 lim, c/dope = 39.5 fim, d m f = 30 
fim (A e ff = 707 fim 2 ), d c ] a i = 329 fim. The fiber was coiled with diameter of 53 cm. The 
signal seed was 30 W, the pump power was <1000 W at 976 nm, in a counter-pumped 
geometry. The measured threshold was at an output signal power of 380-500 W. 

Eq. (1) suggests a gain peak at 2420 Hz, but hole burning would reduce this to nearer 
the measured 2000 Hz. The agreement between measurement and model is reasonable 
in this case. 

11. Conclusions 

We showed using our detailed numerical model that the frequency of the highest mode 
coupling gain behaves in expected ways. If we ignore transverse hole burning so the 
oscillating portion of the heat profile is proportional to the product of fields in the two 
coupled modes, the predicted frequency scales as 1/A e g or l/^ ff as anticipated for a 
diffusive thermal wave. Furthermore, Fm is quite close to the value given by Eq. (3) 
which is an estimate based on thermal diffusion across a heated cylinder of radius r e g . 



We also showed that transverse hole burning tends to move the distribution of os- 
cillating heat away from the axis of the core which tends to decrease Fm- Other effects 
tend to increase Fm- For example, confining the doping or reducing the mode size by 
thermal lensing both tend to increase the frequency of the gain peak. Coiling the fiber 
or cooling only one side both push the modes toward the outside of the bend and away 
from the cooled side, which tends to compress the modes and increase Fm- 

Unfortunately, comparisons between measurement and model are not yet sufficiently 
advanced to verify these predicted effects. The measurements in most cases bear signs 
of input pump or signal modulation, and this means that the modulation spectra may 
strongly influence the measured oscillation frequencies. The obvious next step is more 
careful measurement of seed and pump modulation spectra. 



